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Abstract: Aiming at a model-independent analysis of possible new physics effects in 
semileptonic processes at various energy scales, we list and study a complete set of SU (3)c x 
SU{2)l X U{1)y invariant 4-Fermi operators which consist of a pair of quarks and a pair of 
leptons above the electroweak symmetry breaking. We give a full 1-loop renormalization 
group treatment of the evolution of the Wilson coefficients associated with these 4-Fermi 
operators between low energy (~ meson masses) and the cutoff scale A, ~ (1 — 10) TeV, 
where we assume new degree of freedom beyond standard model will begin to appear and 
an ultra-violet completion of our effective theory will take place. 

Motivated by the existing phenomenological bounds, we argue that the new CP violation 
can only stem from the scalar and tensor types of 4-Fermi interaction. Some interesting 
constraints are obtained by studying the universality of kaon and pion leptonic decays, CP 
violating polarization of K'^^, charged lepton anomalous magnetic moments, and —> ey) 
like rare decays. In particular, we can use the limit of electron dipole moment to constrain 
the size of the CP violating triplet correlation in the e'’'e“ —> tt process. 
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1. Introduction 

Intense experimental efforts in the search for new physics have resulted in establishing 
the Standard Model (SM) as a very good effective theory at the weak scale given by the 
Higgs boson vacuum expectation value of u ~ 250 CeV and below. Yet there are strong 
arguments that it is one with a cutoff scale A much lower than the Planck scale; perhaps 
as low as a few TeV. More importantly the observation of neutrino mixings have given the 
first hint new physics beyond v. The exact form of physics beyond the SM is unknown. One 
can proceed by making a guess at this new physics and engage in constructing consistent 
models. The other approach is to make use of the effective field theory (EFT). The main 
assumption here is that at energies below A physical observables are largely insensitive to 
the unknown new physics. The effective Lagrangian is then a sum of the SM term and 
non-renormalizable ones which are the results of integrating out the unknown degrees of 
freedom. This bottom up approach clearly has its drawbacks. One of them being that 
there can be many such non-renormalizable terms. This is not as hopeless as it seems. 
This operators can be classified by their dimensions and the higher ones are suppressed by 
higher powers of 1/A. Strictly speaking the cutoff scale for each set of operators need not 
be the same. Hence, we need only focus on the lowest dimension operators. 

The EFT approach is built upon two crucial ingredients. The first one to correctly 
identify the symmetry operative below the scale A. The second one is to know the degrees 
of freedom. The phenomenal success of the SM in confronting experiments suggest strongly 
that we take the gauge symmetry of the SM to be operative between A and v. Below the 
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electroweak scale the gauge symmetry is SU{3) x 17(1). To this we may also incorporate 
baryon number so as to ensure proton stability. This is not mandatory as is well known 
that operators leading to proton decay can be suppressed by mass scale of order 10^® GeV. 
While taking the SM as the gauge symmetry below A is largely not controversial the same 
cannot be said about what states to include in constructing an EFT below A. The most 
conservative paradigm is to take only the SM fields of 45 chiral fermions plus the gauge 
bosons and one Higgs doublet. The effective Lagrangian will then take the form 

£eff ^ _|_ —Ts + -^^6 H-• (1-1) 

Here £5 is the dimension 5 operator constructed from the neutrino and Higgs fields which is 
responsible for generating Majorana neutrino masses for the active neutrinos. The seesaw 
scale is denoted by A' which may or may not be the same as A depending on the origin 
of neutrino masses. If this is non-vanishing then lepton number is not conserved. On 
the other hand neutrinos may be Dirac particle in which case £5 vanished^. The details 
of how neutrinos get their masses will not be important and we only note that active 
neutrinos are massive and they mix. Ce is a sum of dimension 6 operators composed of 
chiral fermions, gauge bosons, and the Higgs field. The dots are higher dimension operators 
suppressed by higher powers of 1 /A. The number of operators in Cq is over 20 [Q] even 
after the use of equations of motion and not counting family dependence. Clearly a purely 
phenomenological analysis will be unwieldy. It will be a more tractable problem if we select 
a subset of these that are closed under renormalization to 1 -loop and analyze their physical 
effects. 

We begin with the dimension 6 four-Fermi operators that are made up of a pair of 
leptons and a pair of quarks. We are especially interested in the ones give CP violating 
(CPV) effects*. There are 10 such operators which we will divide into two groups: the 
vector type, and the scalar and tensor type. The vector terms are explicitly, 


= (gVQ^)(LS^L^), (1.2) 

= (QaS^Q6')(4S/.ia") , (1.3) 

= (QVQ^)(eS^eO, (1.4) 

(1.5) 

= (nVn^)(LS^L^), (1.6) 

= (JVd^)(eS^e^), (1.7) 

= (uW)(eSeO, (1.8) 

and the remaining scalar and tensor terms are 

Osf' = (1.9) 


Tt is sufficient to have a U{1)b-l symmetry to ensure proton stability and no Majorana neutrino mass. 
However,we are mindful that Planck scale physics is likely to break such global symmetry. 

‘Since we are considering only semileptonic decays and lepton dipole moments the 4-quark operators 
will not be included here. 
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( 1 . 10 ) 

( 1 . 11 ) 


O5?' = {QaU^){Lb’^e^)e'^\ 

= {QaCr^''u^){L}!'cr^^e'‘)e°‘^, 

where Q, L, e, d, u are respectively the left-handed quark doublets, left-handed lepton dou¬ 
blets, right-handed charged leptons, down and up-type quarks of the SM. a, b are SU (2) 
indices and are family indices. We have also suppressed the chiral projection oper¬ 

ators. We use the convention that repeated indices are summed over and our metric 

is (H-) with { 7 ^, 7 ^} = . The purely leptonic 4-Fermi operators are recently 

studied in detail in and we shall follow the notations given there. Then in the weak 
basis we write 

7 2 

= C'^aO^a + C^sfofA + + h-c. ( 1 . 12 ) 

A=1 A=1 

In general the Wilson coefficients denoted by C’s depend on the family index. Also her- 
miticity requires the 7 vector Wilson coefficients to satisfy 

= (cf^y _ (1.13) 

However, for the scalar and tensor types there is no such requirement. Clearly the most 
general flavor structure is too unwieldy for a fruitful analysis and we shall make phe¬ 
nomenologically motivated simplification that hopefully reflects some dynamics of the yet 
unknown new physics. 

Even after going over to the mass basis of the fermions there will be flavor changing 
neutral currents (FCNC) generated at 1-loop. These are known to be highly suppressed. 
Following the discussions in Q we make the ansatz 

where are constants that depend on the quark and lepton families. The operator 

Ov 2 contains charged weak currents. Universality tests of tau and Pion decays suggest that 
their Wilson coefficients should be independent of lepton flavors. Thus it is reasonable to 
assume 

C^vf = C^V2- (1-15) 

We shall see in the next section that renormalization will mix Oyi and 0^2 and hence 
to satisfy the above mentioned experimental constraints we take Cyf^ = Cy^. Thus, we 
are emboldened to extend it to all vector Wilson coefficients and drop the dependence on 
the families indices. In fact most four dimensional unified models have this feature. Some 
notable exceptions are provided by extra dimension models and split fermion models |^. 
With that in mind we rewrite the effective Lagrangian in the mass eigenbasis as 

7 2 

-Ce = J2 CvAO^f + Y, C^aO^a + + h.c. (1.16) 

A=1 A=1 

From now on all our operators expressed will be in the mass eigenbasis. Notice that 
the flavor insensitivity of the vector operators are not extended to the scalar and tensor 
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operators. Certainly the Cs^t^s in Eq.( 1.16 ) differ from those in Eq.( |1.12] ) by products of 
rotational matrices that diagonalize the fermion masses. For simplicity we have kept the 
same notation. 

There are eight dimension-6 dipole operators with two fermions. They are listed as 


DY = , (1.17) 

(1.18) 

, (1.19) 

dH = ( 1 . 20 ) 

, ( 1 . 21 ) 

Di^ = HQ^a^^'^d^B^, , ( 1 . 22 ) 

(1.23) 

(1.24) 


where H = and B, W, G are the field strength of U{1)y, SU{2)l and SU{3) respec¬ 

tively. 

Notice that the operators are in general not flavor diagonal. After electroweak symme¬ 
try breaking (EWSB) we can substitute H (v + h^)/\/2. The terms from Eq.(1.17, 1.18|) 
become 

^(eV^"i?e^F^,) + h.c. (1.25) 

V2 


rD 


They will induce flavor violating or conserving radiative transitions that result in lepton 
electric dipole moments (EDM’s) , {g — 2)^, /i —> ey, and 6 —> sy at the level of the Born 
term. All these are strongly constrained by null experimental observation beyond SM |^. 
For example, the muon g — 2 constrains Tie and electron EDM limits Tm and they 
have to be < 10“®. ^From this observation we shall adopt a phenomenologically motivated 
assumption that any reasonable new physics beyond SM will only induce negligible tree- 
level dipole operators. 

Proceeding further one recognizes that the dipole operators can be induced at 1-loop by 
the 4-Fermi operators and SM interactions. In Figjl| we show the possible 1-loop diagrams 
that will mix the 4-Fermi operators and the dipole operators. The first three terms are 
vanishing from the above discussion. Since only the t-quark enjoys 0(1) Yukawa coupling 
and thus the mixing involving other light quarks are severely suppressed by their small 
Yukawa couplings and can be neglected. Fig.2(d) shows that we only need to consider the 
dipole term induced from the tensor operator involving the third generation quarks, Oj, ’ . 
Thus, we only need to consider the first two dipole operators involving 2 leptons. At low 
energies for EDM studies the relevant operator involves the photon thus we define 


+ h.c. 


(1.26) 


and the Wilson coefficient is = c^^^cos0^ + in standard notations. The 
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Figure 1: The 1-loop diagrams which mix up the dipole and 4-Fermi operators. Some representa¬ 
tive chirality configuration are shown. AjB stand for the SU{2)l/U{\)y gauge boson. 


complete effective Lagrangian 

7 2 


ai Liie wears. 




A = E CvAO^f + Cg^og^ + h.c. 

A=1 A=1 A=1 


(1.27) 


The physics is clearer if we expand Cq in terms of their components. We split it into 
the charged current (CC) terms and the neutral current (NC) terms. After some algebra 
the NC terms are given by, (L = ^-^,R = 


+ ^7"^ 


-h 


{Cvi + Cv 2 )L + +Cy5^] u' X 
CviL + CvaR] X 
CviL + Cys^] X (eS^Ae^) 
(Cyi + Cv 2 )L + CvaF^ d} X 


+ ^7'^ 

-h M*7'^ 

+ &Jf(XRd^){e^Le^) + &Jf{lfRu^){e^Re^ 

+ Ru^){e^aij_yRe^) + h.c. 


Cv^L + CyjR 
CvsL + CvgR] d* X (eS/ri2e^) 


rt* X (e^7^i?e^ 
f—k_. f>^k\ 


(1.28) 


If one probes the effective neutral current couplings using neutrinos versus electrons 
one would find the corrections to the SM are very different. Interestingly the NC effects of 
the electron can contain terms with different Lorentz structure then that of the SM. On the 
other hand the new NC couplings of the neutrinos are only of the vector and axial vector 
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types. This follows from the assumption that there are no right-handed neutrinos below 
A and the gauge symmetry of the SM. The structure of Eq.( 1.2S| ) also points to polarized 
electrons being powerful and versatile probes of new physics. 

The CC terms are explicitly: 


- = c'y2{ddLd){d + &^{^\dRd){dLJ) 

-&Jf(XRd){v’^Re^) - Rd){v’^af,^Re^) 

+h.c. 


(1.29) 


Cy 2 differs from its NC counterpart by quark and lepton rotational matrices. In principle 
if the parameters involved are all measured in the future the relation between Cv 2 and Cy 2 
will be a test the above framework. Currently we only have limits on these parameters. 


The Eq.(1.2£) renormalizes the canonical 4-Eermi interaction; thus Cy 2 can be revealed 
by precision measurements comparing tau and muon leptonic decays with neutron beta 
decays as done in classic universality tests. The scalar and tensor terms adds incoherently 
and hence are not important. Nonetheless, these latter terms are very interesting since 
their discovery will be strong indication of new physics. They can be probed by K and 
B meson decays; for example by measuring the energy spectral of the charged leptons in 
semileptonic decays. This is analogous to what we have studied previously in /r decays 
Furthermore, the Wilson coefficients are in general complex and thus can give rise to novel 
CP violation effects via interference with a SM amplitude. A prime example will be muon 
polarization measurements in which are unique low energy probes of these 

terms |Q]. Details will be given later. 

Now we return to more theoretical issues. The Wilson coefficients are defined at the 
scale A. In a top down approach one would be able to calculate them in terms of the 
parameters of the new physics theory. The bottom up path taken here does not afford 
such luxury. We can only treat them as parameters to be determined or constrained by 
experiments. To date all relevant CPV experiments are done at energies much lower than 
A. We must then evolve the Wilson coefficients from A to u and then to the scale of K, B 
mesons or charged lepton masses. This is done by solving the renormalization group (RG) 
equations for the Wilson coefficients and matching the boundary conditions as one crosses 
each mass scale. The results here can be used to relate CP violation quantities measured 
at high energies at colliders to those at lower energies and vice versa. As an example some 
leptonic processes cases are studied in |^. 

In this paper we are more interested in processes below the weak scale. In the next 
section we give a discussion of the evolution of the Wilson coefficients from A to the weak 
scale and calculated the anomalous dimension of operators. Below the weak scale we 
integrate out the t-quark and the heavy SM gauge bosons. The RG running then will 
involve only U{l)em and SU{3)c and will be discussed in detail. We choose the low energy 
boundary condition at rup = 1 GeV. In section III we use the effective Hamiltonian and fold 
in the appropriate form factors to estimate the strength of signatures for CPV experiments 
in K meson and pion as well as semileptonic decays, and the electric dipole moments of 
lepton and neutron. They give stringent limits on the scalar and tensor coefficients. As an 
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example of the synergy between precision low energy measurements and high energy collider 
experiments we calculated the CPV triple momentum-spin-spin correlation in e'’'e“ t + i 
at the linear collider. Our results are summarized in the concluding Sec.IV. 
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Figure 2: The connection between high and low energy interactions in effective field theory. The 
symmetry at each mass scale is given. All notations are given in the text. 

A schematic representation of this approach is given in Fig.(^). Although we have 
made the very conservative assumption of no new states below A the effective operative 
approach can be easily generalized. For example if there were sterile neutrinos below A the 
list of additional lepton operators is given in |^ . If the minimal supersymmetric standard 
model turns out to be the correct description of physics just beyond the standard model 
then A can be identified with the supersymmetry breaking scale. Moreover, the operator 
list will have to be extended to include supersymmetric state or states and care has to be 
taken to respect the i?-parity. 


2. Effective Hamiltonian for Leptonic and Semileptonic Meson Decays 


We begin by calculating the quark level effective Hamiltonian for di ujlv. This given by 
the low energy SM with the W boson, Higgs boson, and heavy quarks integrated out plus 
the relevant terms in Eq.( p..29 ) with the Wilson coefficients taken at some energy // below 
the weak scale. The term C'y is expected to make a undetectable correction to the Fermi 
coupling constant Gp since we are taking A to be TeV or higher. Thus, to leading order 
in weak interactions and the new physics scale A the effective Hamiltonian is given by 


= ^^A{n){VijPikua^^Ldj){laf,Luk) 

H—^^^2— -{uiRdj){lkLvi) - *^^2 ^ {diRuj){vkRk) 

- '^A2 ^ {dicr^''Ruj){i?kCrf,uRh) 

+h.c. 


( 2 . 1 ) 














where Vij is an element of the CKM matrix and Pik is an element of the neutrino mixing 
matrix. The Wilson coefficient A is unity in the free quark limit and are modified by SM 
gauge interactions as ones goes from A to the electroweak scale and then further down. 
The QED and QCD running of A(/i) is well known whereas similar behaviors of scalar and 
tensor coefficients, to the best of our knowledge, have not been presented before. Hence we 
shall briefly discuss below how one can obtain the leading logarithmic (LL) contributions 
by renormalization group methods. This also serves to establish our notations. 

The operators Osy,T,D which we shall generically denote by 0\ are bare operators. 
Upon renormalization they will mix via the equation 

10 

oi = ^ z-/^{ZLT/\Zeri\ZQri\z,ri^ot (2.2) 

B=1 


where Zg, Zq and Zq {q = u, d) are the wave function renormalization constants for the 
various fermion fields, n, m, r, s are the number of such fields in each of the Oa operators. 
Thus, n, m,r,s = 0,1, or 2. The sum here runs over the 10 terms of Os, Oy and Ot and 
the prime denotes renormalized quantity. 

The renormalized operator will depend on the t’Hooft renormalization scale fj, 
whereas the bare operators do not. Correspondingly the ^-dependence of the Wilson coef¬ 
ficients will be such as to render the renormalized effective Lagrangian independent of /r. 
This leads to the renormalization group equation (RGE) for the coefficients in Eqs.( Ll^ ): 


H-yCA + IAbCb = 0 


(2.3) 


where '^ab is the anomalous dimension matrix which is non-diagonal. The values of Ca at 
any scale m are obtained by solving the above equation plus boundary conditions which 
are the values given at another scale say the weak scale M\y. At present there are only 
limits on a few of these coefficients from LEPII measurements. We now return to discuss 
'yAB calculation. 

Below the weak scale the renormalization of the operators are given by photon and 
gluon exchanges. We will ignore all Yukawa couplings except that of the t-quark. The 
calculation becomes the evaluation of the high energy parts of the depicted diagrams in 
Fig. i The wave function renormalization graphs also must be calculated but not shown. A 


(a) 







Figure 3: 1-loop corrections to the 4-Fermi effective operators by SM gauge bosons represented 
by the wavy lines. 

lengthy calculation gives the following result for the 4-Fermi operators anomalous dimension 
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matrix: 


/Gi 




\ 


2ai 

0 

0 

0 

0 

0 


0 0 
—2ai 0 

0 4a 
0 
0 
0 


0 

0 

0 


0 

0 

0 

4q!i 

0 

0 


0 

0 

0 

0 

—8ai 
0 


0 

0 

0 

0 

0 

Sas 


G 2 J 


in the basis of {cyi,,y7, C51, CS2, ct}- And Gi ,2 are the two by two matrices: 


(2.4) 


Gi 


G 2 


—ai — 302 802 

602 —oi — 7 o 2 


( —^oi — 803 SOoi — 602 
|oi + §02 |ai + §03 


(2.5) 

( 2 . 6 ) 


After diagonalization and putting this into Eq. ( ^.31) and solving it gives the running of the 
scalar and tensor Wilson terms between A and M^y. 

To get a qualitative feeling of the RG running above the electroweak scale we first take 
the running of SU{2) and U{1) to be much smaller than from SU{3) when the cutoff is 
below 10 TeV. This is in accordance to usual expectation. Considering the SU (3) running 
alone leads to following simple solution. Suppressing indices we define the enhancing factor 




G(A) 

G{Mw) 


(2.7) 


Then we have 

^yi..V"7(A7vV) A) ~ 1.0 ± 0(0.01) (2-8) 


since only U{1) and SU{2) running play a role here. On the other hand QCD running 
affects the other coefficients and we have 


^5i(A7iv, A) ~ Qs2{Mw,^) ~ 




/ a3{Mw) \ 

V «3(Mt) J 


K asiMt) J 

a3{Mt) \^ 

as (A) J 


11 

23 


/ asjMt) 

V as (A) 


_ 4 _ 

11 


(2.9) 

( 2 . 10 ) 


At A = 1 TeV, we find that ^ 51,2 = 0.865 and Gt = 1-049. We can also give approximate 
expressions for the above quantities : 


GsI,2{Mw-,^) 


0.0456 (!^) 


0.036 / A \ 

, Qj'{My/, A.) ~ 1.0188 ( — ) 

\mtj 


0.012 


( 2 . 11 ) 


where the number factors in front of the parentheses are the result of RG running from 
My/ to mt = 174 GeV. 

Below the EWSB scale the operating gauge symmetry is [/(l)em x SU{3)c- It is well 
known that there are no large LL contribution from QCD for semileptonic decays in the 
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SM. The perturbative QCD corrections have been calculated and is found to be small 
henceforth, we shall neglect it in A{^). However, QED gives rise to LL corrections in 
H(/i) 1^ and the beta function is obtained from a subset of the diagrams listed in Fig.|3[ 
Explicitly we have 


A{mi,mw) 


a{mc)' 

3 

8 

a{mr) 

9 

32 

aimb)' 

9 

38 


aimi) 


airric) 


a(mr) 


a(mb) 


( 2 . 12 ) 


where i = d,s. In the above we have integrated out simultaneously the t-quark and the 
gauge bosons and then the light fermions in succession, is either the strange quark 
mass or the proton mass and we have noramlized A(mi,mi) to unity. Using the values 


a-^{Mw) = 128.0, 
a~^{mb) = 132.14, 
a~^{mr) = 133.33, 
a~^{mc) = 133.69, 
a~^{mp) = 133.91, 


(2.13) 


where we choose mb = 4.3 GeV, me = 1.3 GeV, and mp ~ 1 GeV. We find 

A{mp,Mw) = 1-0107, (2.14) 


which is an important correction for precision measurements. 
Eor Csi, we hnd 

i (r “ ’ 

and have the solution: 


Csiirrip) 


a{mc) 

1 

8 

a{mr) 

3 

32 

a{mb) 

3 

38 

a{Mw) 

_a{mp) 


a{mc) _ 


a(mr) 


a(mb) _ 


X 


'asime)' 

9 

asimbY 

25 

az{Mw) 

asimp) 


asimY) 


a^{mb) _ 


12 

23 


Csi{Mw) ■ 


(2.15) 


(2.16) 


Using the boundary condition as{Mw) = 0.120 and the standard 1-loop QGD beta func¬ 
tion, we get 


asimb) = 0.211, 


as{mc) = 0.316, 


Oi^Ynrip) = 0.359, 

(2.17) 

Csiirrip) = (0.9965) x (1.7232) x CsiiMw ), 

(2.18) 


where the first and second brackets are the QED and QCD RG running effect respectively. 
We see that the RG effects enhance the coefficient almost by factor 2. 


- 10 



































The running of Cs 2 and Ct are more complicated since they are coupled: 


a 


4/3 8 


dvr \ 1/6 -40/9 



(2.19) 


in the basis of {Cs 2 : 

An analytic solution for Eq.( p.l9D is unavailable. For a solution numerical methods 
are required. However, since the QED effect is expected to be much smaller than QCD; for 
simplicity, we can just keep the QCD part to have a rough idea how the Wilson coefficient 
evolves. Using the same parameter setting as in Eq.(2.18|), we get 


Cs2{mp) ~ 1.72Cs2{Mw) , CT{mp) ~ 0.83Ct{Mw) 


( 2 . 20 ) 


which should be sufficient for order of magnitude estimates of new physics effects. It is 
interesting to note that the scalar coefficients increases as the energy decreases whereas the 
tensor one does the opposite. 


3. Phenomenology 

3.1 Leptonic Decays of K and vr Mesons 

We can now apply the above formalism to leptonic decays of pseudoscalar mesons. In 
particular the branching ratios Rt^ = B(tt —> ei^/vr —> ixv) and Rk = B{K ev jK 
which are precisely measured. Theoretically they have been precisely calculated in the SM 
§ with the uncertainties below the experimental limits; thus making these decays very 
sensitive tests of new physics and in particular scalar interactions. 

To calculate the amplitude a pseudoscalar meson M (7r+, A'+oril+) decays into a 
lepton pair we need to introduce the following form factors 

{Q\ua^-i5dj\M{p)) = ifMP^ , 

{0\ua5dj\M{p)) = ifM /{mi + rrij)) , (3.1) 

where M is the meson mass. The pseudoscalar form factors take the values fn = 130.7(4) 
MeV and fx = 160(2) MeV although we do not need these values. On the other hand fs 
is not known but expected to be about 200 MeV. 

Since the pseudoscalar meson carries only one momentum, p. It can’t have an anti¬ 
symmetric form factor due to Lorentz invariance. 


{0\uia^’'dj\M{p)) = 0 , 

(0|u,u^"75dj|M(p)) = 0. (3.2) 


However, one can have nonzero form factors for tensor interaction if there is an extra 
particle in the final state such in the vr” —> ei"y case |1T| to which we shall return at the 
end of the section. 
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Substituting Eq.(^) into Eq.(^) the amplitude for the decay dj —> Uilu^ for a specific 
charged lepton I but sum over all three active neutrinos is then given by 






VijA{mi) i^Piklp^'yij.Luk 

V2M^ 


- C''d^^*{mi)) ( ILu, 




4GirA2(mi + TTij) V 
and notice that only a sum over k is taken. From this we obtain 

KI-K!: 


(3.3) 


R^ = Rt^ X [I + 


GfA^ 


where 


jSM _ 




a 


ml{M^ - mlf 
is the SM prediction Q, and 

1 m2 


1 - 16.1- = 1.2354(2) X 10 

TT 


-A 


KL = 


y/2A{mi, Mw)Vud m{mu + md) 


[Pki ( 


^kl[^S2 ~^S1 


(3.4) 


(3.5) 


(3.6) 


Because we have factored out the SM helicity suppression factor in Eq.(3.4) the correction 
factor given by Eq.( |3.6D has a Ijmi. This is a reflection of the helicity flipping nature of 
scalar interactions. 

Compared with the data R^^ = 1.230(4) x 10“^|^ which is consistent with the SM 
we have the limit: 

ly-e _ 

" < 0.007 (68% C.L.) (3.7) 


GpA^ 

Besides the overall scale factor A for the new physics, Eq. (|3.7D gives valuable informa¬ 
tion of these scalar coefficients. The first observation is that the phases of G^^^^ and G^^^^ 
must not accidentally cancel for it to be useful. Next we define 


^7T,kl ^11,kl* ^11,kl /I \ 


(3.8) 


and consider the following general possibilities which cover most of the relevant models for 
this problem: 


A. There is no hierarchy (NH) between the electron and muon modes. 

An an example we take 

Then dominates in Eq.(|3.7|) and taking MT^/{mu + md) = 15, we have: 


E 


IHe 


ke 


G 


11,fee* 
S 2 


^ll,ek\ 
^S1 ) 


< 2.8 X 10“^ 

NH 



(68% C.L.) 


(3.10) 
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B. An exact hierarchy exists between the electron and muon modes. 

As an illustration we take the coefficients to be proportional to the electron and muon 
Yukawa couplings, ye ~ 3 x 10“® and ~ 6 x 10“^ respectively. Such is the case for 
two Higgs doublet model (2HMD). Explicitly we have 

cf^ = yeC, 

, (3.11) 


where C is an undetermined constant which for simplicity we take to have no depen¬ 
dence on the neutrino indices. In this case the constraint from Eq.(^) becomes 


[Pke 

k 




(3.12) 


Clearly the constraint from lepton universality tests is not as stringent as the previous 
case; nevertheless is still an interesting one. Furthermore, it also demonstrates the 
importance of taking into account neutrino mixings. 


The above analysis can be repeated for K decays. We have 

X (l + 


GfA2 


where is the SM prediction and 


K^k = 


Ml 


Mw Wus ^t{^u + rria 




^kl [ ^52 “ ^51 


The most recent data from NA48/2 |]T] is 

= (2.416 ± 0.043,tet ± 0.024,j,,t) x lO”® 


(3.13) 


(3.14) 


(3.15) 


which improves upon the PDG value = 2.44(11) x 10“^. This is to be compared 

with the SM value of = 2.472(1) x 10“^ |^]. Using Mif/(m„-|-ms) = 3 and 14s = 0.224, 
we have the limit: 




^21,fee* ^12,ek 

(_x QQ (_y 


ke I ^S2 


'SI 


< 6.4 X 10 


-5 


NH 


For the exact hierarchy case we get 

[Pke-Pkf.] 


< 21.3 




TeV/ 


(3.16) 


(3.17) 


where C replaces C of Eg.( 3.11 ). 

As remarked earlier we have to use the radiative decays to get constraint on the tensor 
coefficients. However, as has been pointed out by due to QED correction the tensor 
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operator will induce scalar operator. And it’s been widely believed that the limit on Ct 
from 7r“ —> ev is two orders better than the direct constraint from vr” —> ev^. 

Since we have a RG improved result given in Eq.( p.l9 ), it’s interesting to see how this 
affects the above conclusion. To that end we carried out a full numerical analysis for the 
coupled equations Eq.( 2.19|) and the result is displayed in Fig.^ where the figure caption 
explains our notations. 



Figure 4: The full RG running for the coupled scalar and tensor Wilson coefficients from Mw to 
rup ^ 1 GeV. The x-axis is the ratio r = CsiICt at Mw, and the left-handed y-axis is the ratio 
r at TOp. The dash line (blue) is the amplitude at low energy, we set A = v^|Cs 2 p + |GtP = 1 at 
Mw, refer to the axis at the right-handed side. 


Assuming there is no hierarchy between the Wilson coefficients Cs 2 and Ct at EW 
scale, for instance the two are within one order of magnitude. Basically, we find the RG 
evolution can be simply summarized by Eq.( p.2Ci ). 

Let us consider the extreme examples that {Cs 2 ,Ct) = (0,1) and (1,0) at EW scale. 
The first corresponds to only tensor interactions are induced by new physics whereas the 
second has only scalar being generated. Our numerical solution gives (—0.03, 0.84) and 
(1.717, —5 X 10“^) at rrip respectively. In other words, if the new physics beyond SM gives 
only tensor interaction at the weak scale; then at meson mass scale the RG evolution gives 


Cs2{mp) 

CT{mp) 


-0.036. 


This is to be compared with the estimation given in |12] 


Cs2imp) 

Crimp) 


--In^ 
TT mi 


- 0 . 020 . 


(3.18) 


(3.19) 
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We see the full RG running gives about 80% correction. Using Eq.( 3.10|) and Eg.( 3.18 ) we 
get the upper limit 


|ineCrK)|< 7.7x10-^ 

which is a factor of ~ 250 better than direct measurement of vr ^ ewj. 


(3.20) 




3.2 Polarization in K~^ 

It is well known that the transverse muon polarization in semileptonic pseudoscalar mesons 
decays are sensitive CP violation tests of effective scalar interactions |13, 14|. We revisit 
these studies by adding the RG running of the effective operators to the previous studies 
which were conducted mostly in the context of specific models II HI- 

The effective Hamiltonian at the quark level is the same as in Eq.(p.l|). The form 
factors for K^{pk) —> decay that we need are 


(7r°|s7^u|R:+) = f+pl + f.p't 
{Tr^\su\K~^) 


- m2 


m.. 


• /r 


(7r°|sa^‘^n|R:+) ~ i— , 


(3.21) 


where p± = {px ^P-w) and we have neglected the u-quark mass. Again, the form factors of 
(s 7 ^ 75 u), (s 75 u), and (sa^^j^u) vanish due to parity. The form factors are functions ofp?_ 
and is normalized to /+(0) = 1. Furthermore, one usually defines ^ In the second of 

Eq.( |3.2lD we have omitted a term proportional to ^ by taking advantage of its small value; 
i.e. 1^1 = 0.124 1^]. In the third of Eq.( |3.21] ), we have dropped a term Sfr x (e^^^^Pi^oPw/s) 
which can be taken into account by redefining fx (/t ~ < 5 /r) when contract with the 
lepton current part. The decay amplitude is 

-M= V2GFf+ \ V:,P2kA{ms){p'^ + CP-)i>7aLp 


+ 


{Ml - Ml) 
2^/2GFJM‘ms 


— i 


r ^o21,fc2 / \ 

Jt Cj.' {ms) 


/+ \/2GfA'^Mk 


- G^S2 i^s)] 
PxPn^kCTapRp ■ 


(3.22) 


The above can be easily compared with the standard notation y]. We further define two 
variables 


^5 = 


1 


Ml-ml Gjf 


= - 


2^/2VusA{ms) Mxms IGGf ’ 

h/U cl 


where 


VWusA{ms) A^Gf ’ 


(3.23) 

(3.24) 

(3.25) 
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and 


(3.26) 




The CP violating transverse muon polarization is found to be: 

o2 


Pl 


^ . PK-{Pu- Pfi) + ml/2 

-’"’f* +-Mj- 




Pli X pi, 


(3.27) 


where jPs are leptons’ 3-vector momentum in the kaon rest frame and we have ignored the 
correction of higher power in and 


= {2E^Mk - ml) ( 1 - 


E-k + Efj 

Mk 


1 


- + m^-m - 2 E^Mk) 1 - 


(' 


Plugging in the numbers, we have 
Pi 


Pk ■ (Pu — Pu) + mf,/2 „ 

0.383„,Cf - 


_mp_ 

AMj. 


TeV Y p^xpy 
A ) $ 


(3.28) 


(3.29) 


It is interesting to note the imaginary part of both tensor and scalar coefficients are sensitive 
to where measurements are made in the Dalitz plot distribution |17|. Ct has maximum 
sensitivity when — Ei, is largest in the kaon rest frame . From the current limit \Pt\ < 


0.0050 and assuming that /t — f+, we have 

lUmCg^l and lUmC^I < 2 x 10“^ . 


(3.30) 


Notice here that the results from K 12 decays limits the real part of the scalar coefficients. 


3.3 Constraints from Electron and Neutron EDM’s 

It is well known that the EDM of the electron is a very sensitive test of new physics phases. 
This is due to the fact that the SM gives an undetectably small value \de\ < 10“^® e-cm. 
In addition there are several ongoing experimental efforts to reach the level of 10“^*^. If 
discovered it will be a clear sign of new physics. Similarly the neutron EDM, is also 
very sensitive to new physics phases and are usually connected with electron EDM in model 
dependent ways. Hence, they are very complementary tests of new physics. Although the 
SM prediction for dn is not as robust as dg but is generally taken to be < 10“^^ e-cm. This 
is still very much below the experimental bound of 

|d“P| < 6.3 X I0"2® e-cm. (3.31) 

Among the terms in the effective Lagrangian new phases reside in the scalar and tensor 
Wilson coefficients in and Some NC terms will contribute to di at the 1-loop 

level and the CC term will induce it at the 2-loop level. This is explored in detail below. 
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Figure 5: 1-loop contribution to the EDM of a charged fermion / from semileptonic 4-Fermi 
operators denoted by the box and / represents I or a u-quark and f a u-quark or I respectively. 


3.3.1 EDM at One Loop 

The generic one loop diagram for dj is given in Fig.(|^. It has been shown that the RG 
running of di below the EW scale is not significant and we shall ignore it. For the quark 
EDM, QCD correction is the biggest contribution. We found its QCD anomalous dimension 
is same as the tensor operator which we have already found its solution in Eq.(2.20|), i.e. 


du{mp) ~ 0.83du{Mw) 


(3.32) 


Hence, for our purpose we can neglect the RG running in both cases. 

It has been shown that scalar operators will not contribute at this level. The vector 
terms in Eq.( [1.28 ) have real coefficients and hence will play no role. This leaves the terms 
and only the u-type quarks come into play. A straightforward calculation leads to 


dp_ — 


eNr, 


3yr2A2 

i=u,c,t 


In 


1717 


(3.33) 


at the EW scale which is where we set the boundary condition for the RGE of EDM 
operator and W is the number of colors. For simplicity, we ignore the top threshold effects 
and the RG evolution between mt and Mp/. Similarly, we have 


du — 


e 

27r2A2 


m; In ^ . 

TTi] 


(3.34) 


Notice that there is no dd at this level due to the SU{2)l symmetry. 

Clearly the most important term is due to the t-quark. This has the unusual behavior 
of depending on the internal fermion mass and independent of nie. The chirality changing 
nature of tensor interactions is responsible. Thus, one expects the same formula to hold 
for dp and d^- with obvious change to the lepton indices. ^From the experimental limit of 
\de\ < 1-7 X 10“^'^ e-cm 0 we deduce that 


JmC, 


33,ee 


< 1.4 X 10“® 



(3.35) 


which is much more stringent than one obtains from meson decays; see Eq.( |3.30 ), albeit 
for different flavor indices. 
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Assuming that the leading contribution to dn comes from du and using the quark model 
relation dn ~ — du)/S for estimation, we obtain from Eq.( 3.31 ) a less stringent bound 


JmC, 


11,rr 


< 8.4 X 10 


-6 



(3.36) 


3.3.2 EDM at Two Loops 

Unlike the complex coefficients in can only induce EDM at the 2-loop level. 
Similarly the scalar coefficients of come into play now. The generic diagrams for 
this are displayed in Fig.(|^ Now both the scalar terms Csi ,2 contribute. From the above 




Figure 6: 2-loop contribution to the EDM of a charged fermion / from semileptonic 4-Fermi 
operators denoted by the box. The internal gauge boson line is either a photon or W-boson. 


discussions we can ignore the tensor terms since they are constrained to be small (see also 
the next section). 

An examination of and the structure of the Feynman diagrams shows that only 
the flavor diagonal terms come into play. Since the lepton masses are small they can be 
neglected in the calculation and we only need to keep the quark masses. Also the leading 
contribution comes from internal photon exchange. Then the NC contribution to the EDM 
di of a charged lepton Z = e, p or r is given by 


dz(NC) = 


ea 


dSvr^A^ 
and the functions are 




i=u,c,t 


m- 


, (3.37) 


= 

rdx 

Jo Jo 1 + y(i - y)x 

(3.38) 

W = 

rda:[\y . 

Jo Jo 1 + y(i - y)x 

(3.39) 


These results extend our previous expression in |^] where only partial contribution was 
considered. We have Tt = {10.872,59.552,150.780} and Tb = {8.596,13.200,17.808} for 
A = {1,10,100} TeV. When z S> 1, they can be approximated by 


Tt{z) — {2 — \iiz + \d?z) , 

J^b{z) — — (2 —Inz) . (3.40) 
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There are two helicity flips in the above result. The first one involves the scalar 4- 
Fermi interactions Csi, 2 - The second one is the helicity flip from the quarks in the loop 


due to the chiral nature of the weak interactions and is explicitly displayed in Eq.(3.37). 
As expected, in the limit of massless fermions there is no EDM. It is also clear that only 
the t and b-quarks are important. 

For the contributions from the calculation is similar. Now the internal gauge 
boson exchange comes from the W boson. The result is 


di{CC) 


ea 




IGvr^ sin^ 

O Jo Jo 


mtJmC ,,,’ Gt ( ) +mbJmC^,’ Gh 


'S2 

(2-3y)(2-2/)x 




Gb{z) = -- 


0 Jo + yx) (1 + x) ’ 


(3.41) 

(3.42) 

(3.43) 


where we have kept only the t quark mass in the quark propagators. The analytic expression 
for G is rather complicated and not illuminating. From numerical study we have Gt = 
{1.708,11.450, 31.670} and Gb = {-0.0824, -0.347, -0.6343} for A = {1,10,100} TeV. For 
our purpose it is good enough to use the following approximation when z ^ 1: 


Gt{z) — ^ (13.7 — 6.31 Inz + In^ z) 


Gbiz) 




(3.44) 


Combining the above two contributions, we arrive at the estimation: 


ea 


IGvr^A^ 


mtJxnCg2 




Gt 


sin^ 9,, 


+ mb'JxnC, 


33,11 

SI 


3 ^ sin^ 9„ 


(3.45) 


The physics is now clear. For Csi/s 2 ^ due to the helicity, it can only pick up the bot¬ 
tom/top quark mass in the loop. Barring accidental cancellation between the two scalar 
coefficients terms the upper limit of the imaginary parts can be deduced from electron 
EDM experiment: 


3mC\ 


33,ee 
S2 


< 1.5 X 10 


-6 




< 4.7 X 10 


-4 


A 

IW 

A 


(3.46) 


The results from our one and 2-loop study show that EDM’s are sensitive to the 
imaginary parts of the flavor diagonal scalar coefficients. On the other hand the charged 
lepton polarization experiments directly probe the flavor off diagonal terms although at a 
less sensitive level. 
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3.4 Muon Anomalous Magnetic Moment, and Rare Decays. 

The exact same 1-loop calculation leads to the following effective Lagrangian: 


AC ~ 


27r2A2 


In 


A 2 


mt 


) F, 




+ h.c. 


erUrC^’^^ h? 
47r2A2 m2 



F^u + h.c. 


(3.47) 


The flavor diagonal entries in the first term also contribute to {g — 2)^ and off diagonal 
terms lead to the rare decay {g ^ e + ^). Thus, the branching ratio of the radiative decay 
of a charged lepton L normalized to the decay width of (L ^ + hi) is given by |^: 


B{L 


/ + 7 ) = 48 


2 In^^ 

ami ^ 

ttMI (A2Gf)2 


^ C<33Ti 


2 


+ 



(3.48) 


/,From the 
1.6 X 10“® 


limit B{g 
§, we get 


67 ) < 1.2 X 10-^^ B{t gj) < 6.5 x 10 


-6 


c: 


33 ,fie 


a 


33,Til 


+ 


+ 


a 


33,e/i 


33,/ir 


a 


33 ,lie 


+ 


a 


^33,e/i 

rp 


< 4.0 X 10"^® 


< 6.8 X 10“® 


A 


TeVy 

A 


TeV, 


< 1.7 X 10 


-7 


A 


TeV 


and B{t ^ 67 ) < 


(3.49) 


We see that the lepton flavor changing tensor coefficients are constrained to be very small 
unless A > 100 TeV. 

By the same token we obtain the following modification to the charged lepton anoma¬ 
lous magnetic moment 


, mmi , A2 00 II 

Aai = -^ In —ReCi^^ 


7r2A2 


mt 


/,From the current limits |Aae| < 3.5 x 10 and |Aa^| < 254 x 10 


1-11 


(3.50) 
we deduce that 


33 ,ee 






< 1.1 X 10 


-3 


< 4.1 X 10 


-4 


A 

TeVy 


which are less stringent than Eq.(p.49) but still very strong. 


(3.51) 


3.5 Triple Spin-Momenta Correlations in e+e ^ tt 


We give as an example a physical quantity measurable in the process at high energy which 
can be related to low energy constraints by RGB’s. Consider at the TeV linear collider (LC) 
the process e^{p+)e~{p-, Se) —> t(A:+, st)t{k-), where the momenta and spins are specified. 
One can construct a T-odd quantity involving two independent spins and a momentum ; 
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e.g. p- ■ {st X Se) where j5_ is the unit momentum 3-vectors, and Se is the unit 3-spin vector 
in the electron rest frame. We can relate this to de- 

Since the tensor coefficient has been determined to be around or smaller than 


10 ® from the electron EDM, Eq.( 3.35|) , we need not consider it. Only the effective scalar 
operators will be important. The triple correlation is a result of the interference between 
the SM amplitudes of photon and Z exchange, and the scalar 4-Eermi interaction. To 
match the chirality, either the mass or the spin of electron has to be involved. Clearly, a 
LC with polarized electron beam offers higher possibility of probing CPV part of the scalar 
interaction beyond SM. 

A straightforward calculation yields the CPV amplitude, normalized to the photon 
exchanging amplitude, 


II TO — 

+ 11 + 


mt 


1 + — ] {p- - P+) ■ [St X Sej 


2m? 


{k+ - k-) ■ {stx Se) 


(3.52) 


where s = (p+ + p-)‘^ is the Mandelstam variable, and st is the unit spin vector of the 
t-quark in its rest frame. 

In arriving at the above expression, we have assumed that ;§> Mz and sin^ 9]y ~ 
0.25 for simplicity. Putting in all the RG running we have derived, Eqs.( |2.11 , 2.20| ), we get 
the following relation 




TeV 


f A 

= 0-516 X — 


1.964 


(3.53) 


With the electron EDM constraint, Eq.(3.46), we conclude the relative size of CPV ampli¬ 
tude has to be 

/ /~ \ "^ '\T\ 0-036 

< 2.7 X 10-1 X (4) X (^) . (3.54) 

It is very challenging if not impossible to achieve this kind of precision at the ILC. 


4. Conclusions 

We begin with the assumption that the SM gauge symmetry is valid from the some unknown 
new physics scale A to the EW scale and the matter content of the SM. Using the effective 
field theory approach we studied all the dimension six operators which contain a pair each 
of lepton and quark fields. A previously derived result |^] allows us to cast them in the mass 
eigenstates. By a general assumption that the vector Wilson coefficients are constants we 
can eliminate zero order FCNC. This leaves the most the scalar and tensor operators as the 
most interesting ones that can have FCNC at the level of the Born term. They also provide 
new sources of CP violation in both charge and neutral current channels. Furthermore, 
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in the NC channels there can be flavor changing as well as flavor conserving modes. For 
notational simplicity we generically call the corresponding coefficients Cg and Ct- We 
also argue that stringent experimental limits make the Wilson coefficients of dimension six 
operator of the dipole type negligible and permit us to concentrate on 4-Fermi operators. 

An important aspect of field theory is that these Wilson coefficients evolved with energy 
and this can be calculated by the RG method. Once a coefficient is measured at one scale, 
say the EW scale, its value at higher energy is determined by the anomalous dimensions. 
This robust prediction depends only on the assumed gauge symmetry and the known states 
between the two energies. We have calculated the anomalous dimensions of the complete 
set of semileptonic 4-Fermi operators to 1-loop. The RG running of the associate Wilson 
coefficients are solved between the EW scale and A. 

Below the EW scale the conserving gauge symmetries are U{l)em and color SU{3)c- 
We found that the effect due to U{l)em is not large in agreement with previous studies; 
however the QGD corrections to Cg and Ct are large. The combine effect increases Cg 
by almost a factor of two from the scale of M\y to 1 GeV; whereas Ct is decreased by 
approximately seventeen percent. 

As is well known that pure leptonic meson decays are the most sensitive tests of Cg, we 
refined the previous analysis on the constraints by including RG effects. We also found that 
the large logarithms of the RG method change the result on Ct by almost a factor of two. 
It turns out that the most stringent constraint on all the different Cts arise from 1-loop 
effects they induce in EDM, anomalous magnetic moments, and rare decays of leptons. 

The GPV effects from iJmC '5 do not give rise to EDM’s at 1-loop. At 2-loops they 
constraint the flavor diagonal operators; especially those involving the t-quark. Since these 
are suppressed by an extra loop factor the bounds are not as tight as those for Ct- The 
flavor off-diagonal limes' are most directly constrained by polarization measurements in 
semileptonic meson decays such as in Hence, these are very complementary tests in 
these very general framework of effective operator approach to new physics. Eurthermore, 
given the limits obtained at low energies and the RG effects are not large we do not expect 
any of the Cg coefficients to be measurable at the LHG. On the other hand high energy 
colliders may be more sensitive to the vector terms in This warrants a detail study 

which we shall reserve for future investigations. 

Theoretically the bounds we obtained on Cg and Ct have profound implications for 
model building. They imply that if the new physics scale characterized by A were to be 
below 10 TeV any tree-level origin of these effective operators must have small couplings. 
Although it is not impossible this implies some degree of fine tuning is at work. Well 
known examples are the multiple Higgs models [13| and R-parity violating supersymmetric 
models. An alternative and perhaps more natural one will be that they are dynamically 
suppressed. For example they are generated at 1-loop or higher level in the new physics 
model. Many supersymmetric models have this feature. Another possibility will be that 
the suppression is a result of some symmetry at work above A. Family symmetry is an 
example that comes to mind which will suppress off diagonal terms and leave the diagonal 
ones relatively free. 

We are looking forward to more precision measurements and hope that some of these 
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operators will be measured. Our analysis will be useful to unraveling the new physics at 
work in a general way. 
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